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We  suggest  that  the  electron  Involved  in  an  electron  transfer  reaction  | 

t 

be  represented  as  a  wave packet  expanded  in  a  limited  basis  set  of  j 

appropriate  functions.  Saddle  point  optimization  to  determine  the 
transition  state  involves  moving  the  center- of- gravity  of  the  wavepacket 

t 

together  with  the  atoms  of  the  molecular  framework.  | 

I 
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Introduction 

The  theory  of  electron  transfer  between  molecular  or  atomic  fragments 
in  polar  solution,  in  biological  milieu,  and  at  electrodes  has  developed  at 
the  molecular  level  largely  from  the  point  of  view  of  molecular  orbital 
theory  [1,2].  Thus,  for  example,  one  considers  a  transferable  electron 
initially  located  at  one  atomic  or  molecular  center.  After  successful 
transfer,  the  electron  is  located  at  a  second,  well-defined  center.  In  the 
initial  and  final  state,  the  state  of  the  electron  can  be  described  in  terms 
of  local  state  functions:  xi  and  *f.  The  effect  of  the  solvent  environment 
is  represented  as  a  collective  degree  of  freedom  and  then  separated  from 
molecular  electronic  and  vibrational  degrees  of  freedom  according  to  the 
precepts  of  the  Born-Oppenheimer  approximation.  As  one  of  us  noted  some 
time  ago  [2],  when  carried  out  carefully  according  to  prescription,  the 
Born-Oppenheimer-Holstein  separation  [3,4]  applied  to  the  electron  transfer 
system  yields  a  representation  of  the  underlying  electronic  degrees  of 
freedom  which  closely  corresponds  to  the  familiar  Hartree-Fock- (Roothaan) 

[5]  quantum  mechanics  of  electrons.  In  essence,  the  transfer-electron 
occupies  molecular  orbitals  of  the  reactive  system  in  its  initial  and  final 
states.  This  picture  is  most  clearly  stated  in  Hush's  original  formulation 


[6].  The  state  of  the  electron  can  be  specified  in  terns  of  a  quantity  A 
which  is  used  to  mix  the  initial  and  final  basis  functions: 

♦  -  N(Xl  +  Xxt )  (1) 

and  N  is  the  nomalization.  Adjusting  A  should  fix  the  location  of  maximum 
density  of  electronic  charge  in  the  transition  state.  Indeed,  this  approach 
works  and  is  successful  within  the  context  of  the  electron  transfer  theory 
as  currently  formulated. 

Our  objective,  however,  is  ultimately  to  simulate  electron  transfer 
reactions  in  terms  of  the  discrete  structures  of  solute  and  solvent.  One 
part  of  this  plan  is  to  consider  saddle  point  optimization  of  individual 
configurations  of  particles.  Another  part  of  the  plan  is  to  incorporate  the 
solute  with  its  mobile  electron  in  a  Metropolis/Monte  Carlo  [7]  simulation 
of  the  statistical  transition  state  [8].  Ue  found  when  we  began  this  work 
that  the  conventional  approach  to  the  localization  of  the  maximum  density  of 
migratory  electronic  charge  was  at  best  cumbersome  and  at  worst  impossible 
to  use. 

There  is  a  considerable  body  of  work  on  numerical  optimization  to 
determine  the  transition  state  as  a  saddle  point  [9-14].  In  all  instances 
in  which  these  techniques  are  applied  to  reactions  Involving  the  transfer  or 
rearrangement  of  atoms,  an  assumption  of  the  adiabatic  response  of  the 
electron  is  used.  The  migrating  atom  in  a  reaction  behaves,  for  the 
purposes  of  the  optimization,  as  a  classical  particle;  it  has  a  definite 
extent,  i.e.,  it  is  not  diffuse  as  is  the  electronic  distribution.  Thus, 
the  saddle  point  can  be  located  by  the  manipulation  of  definite  points,  the 
centers-of-mass  of  the  atoms  in  the  system. 

This  approach  to  the  optimization  of  the  location  of  the  electron  in 
the  transition  state  of  an  electron  transfer  reaction  is  not  easily  extended 
with  the  usual  molecular  orbital  view  of  the  states  of  the  electron.  In 


& 


principle,  one  could  optimize  A  of  eq  (1)  together  with  the  positions  of  the 
atoms,  but  this  we  found  not  to  be  the  most  direct  and  transparent  method  to 
use.  Molecular  orbital  methods  are  ideally  suited  to  the  determination  of 
stable,  stationary  states  of  molecular  systems,  states  which  are 
characterized  by  substantial  accumulations  of  electronic  charge  in  the 
vicinities  of  the  atoms  which  make  up  the  molecule.  Thus,  molecular  orbital 
theory  is  not  the  optimum  representation  to  use  to  describe  the 
non- stationary  transition  state. 

On  the  other  hand,  for  the  transition  state  of  the  electron  transfer 
reaction,  we  propose,  first,  that  the  electronic  substructure  of  the 
reactants  be  handled  in  the  usual  manner,  in  terms  of  the  molecular  orbital 
theory.  However,  second,  the  migratory  electron,  should  be  handled  as  a 
wavepacket  which  is  in  turn  expanded  in  a  set  of  basis  functions 
sufficiently  large  to  give  an  accurate  value  of  the  energy  of  the  transition 
state.  The  center-of -gravity  of  the  wavepacket  can  be  manipulated  together 
with  the  locations  of  the  atoms  of  the  substrate  in  order  to  locate  the 
transition  state.  The  wavepacket  is  used  in  the  much  same  spirit  as 
Heller's  [15]  use  of  the  semiclassical  Gaussian  wavepacket  in  scattering. 

Method  and  Example 

We  illustrate  our  proposal  for  determining  the  transition  state  of  the 
electron  transfer  reaction  with  a  simple  model  system.  With  reference  to 
Fig.  1,  consider  two  atomic  centers,  A  and  B,  separated  by  a  distance  R  . 

We  allow  variable  charges  on  the  cores  of  A  and  B,  up  to  two  positive 
charges,  and  a  complement  of  valence  level  electrons,  maximum  of  two.  The 
valence  level  electrons  occupy  simple  molecular  orbitals  of  the  combined  A 
and  B  reactant  system.  Finally,  the  migratory  transfer  electron  is  handled 


as  a  wavepacket,  as  indicated,  expanded  in  a  basis  set  which  consists  of  Is, 
2s,  and  three  2p  Slater  functions  expanded  in  GTO's  [16]. 

The  bo del  system  is  suanarized  in  Fig.  1  as  follows.  The  underlying 
electronic  charge  distribution  which  arises  from  the  interaction  of  the  A 
and  B  species  is  indicated  by  the  light  dotted  lines.  These  lines 
scheaatically  outline  the  regions  of  electronic  charge  accumulation.  The 
transfer  electronic  charge  is  indicated  by  the  bolder  dashed  line.  The 
origin  of  this  distribution,  0  in  the  figure,  lies  in  the  space  of  the  AB 
system.  In  a  saddle  point  calculation,  the  origin  0  would  be  optimized 
along  with  the  surrounding  atoms  to  their  locations  in  the  transition  state. 

tfe  assume  a  closed  shell  or  isolated  valence  electron  system. 

Therefore,  we  ignore  the  exchange  interactions  which  operate  between  the 
migratory  and  the  underlying  electrons.  These  interactions  certainly  may  be 
important  in  many  instances,  but  for  this  initial  report,  we  believe  it  is 
more  important  to  stress  the  major  features  of  the  model. 

The  core  electrons  are  modelled  simply  as  a  single  GTO;  this  yields  an 
electronic  charge  density  about  the  atoms  A  and  B  of  the  type 

P<r)  -  3  *  3  exp(-ra/«a)  (2) 

IT  *6 

In  addition,  of  course,  there  is  an  overlap  charge  density  of  a  similar 
form. 

It  is  worth  noting  that  the  model  system  we  use  can  represent  both 
intra-  and  inter -molecular  transfer.  For  non-bonding,  intermolecular 
transfer,  the  electronic  overlap  distribution  is  vanishingly  small.  The 
height  of  the  barrier  to  transfer  in  this  case  is  essentially  the  difference 
between  the  energy  of  the  free  (or  pseudosolvated)  and  bonded  electron.  On 
the  other  hand,  there  are  examples  of  intramolecular  electron  transfer,  such 
as  the  case  of  the  heptafulvalene  radical  anion  [ 17 ] ,  for  which  the 


interaction  of  the  transfer  electron  with  the  overlap  charge  densities 
should  contribute  importantly  to  the  height  of  the  barrier. 

Results 

« 

Several  calculations  were  carried  out.  The  center -of- gravity  of  the 
transfer  electron  wavepacket  was  moved  along  the  AB-line;  we  calculated 
energies  on  the  AB  line  as  there  was  no  solvent  present  to  warrant 
optimization  off  the  line.  It  is  clear,  however,  that  if  solvent  is 
present,  optimization  will  yield  locations  for  the  center-of-gravity  of  the 
wavepacket  off  the  AB  line. 

The  results  of  our  calculations  are  summarized  in  the  Figs.  2  through 
5.  The  simplicity  of  the  model,  the  absence  of  exchange  interactions, 
precludes  an  antibonding  state  for  the  transfer  electron  for  all 
separations.  Thus,  it  is  possible  to  choose  combinations  of  core  positive 
and  valence  electronic  charge  and  values  for  the  AB  separation  which  yield 
energy  minima  for  the  transfer  electron.  Such  minima  may  persist  even  in 
presence  of  exchange  interactions. 

Fig.  2  illustrates  the  energy  of  the  wavepacket  moving  in  a  fixed  field 
of  two  core  atoms  held  at  a  distance  of  1.2A  apart.  Each  atom  is 
electrically  neutral  at  infinite  separation.  At  this  close  separation,  the 
system  exhibits  energy  minimum  for  the  electron;  however,  the  electron- 
electron  repulsion  at  each  atom  forces  the  transfer  electron  into  the  space 
between  the  atoms  as  being  energetically  the  most  stable.  Note,  there  is  a 
barrier  to  the  penetration  of  the  electron  into  the  atomic  core.  The 
reaction  of  Figs.  2  and  3  corresponds  to  one  of  the  type 

A*  +  A  ^  A  +  A'.  (A) 

Figure  3,  in  contrast  to  Fig.  2,  shows  the  behavior  of  the  transfer 


electron  in  Che  field  of  Che  same  two  aComs,  buc  here  Che  fixed  distance  is 
2.SA.  There  is  an  attraction  of  Che  eleccron  Co  Che  atom- -an  electron 
affinity  Co  form  an  anion,  buc  very  little  penetration  to  the  core.  The 
depression  of  the  energy  at  Che  center  reflects  the  tendency  of  the  energy 
Co  go  to  a  zero  value  as  Che  interaction  with  Che  surroundings  vanishes  with 
large  distance. 

Figs.  4  and  5  represent  reactions  of  the  type 

A+  +  A  ~>  A  +  A+  (B) 

and 

A+  +  B  ~>  A  +  B+  (C) 

respectively.  The  only  distinction  between  these  two  reactions  is  the  fact 
that  for  the  heteronuclear  reaction,  the  parameters  S  and  6  differ  by 
O.lA.  It  is  clear,  however,  in  both  cases  that  there  is  a  barrier  to  the 
migration  of  the  electron  from  one  center  to  the  other. 

Discussion 

We  have  found  only  energy  maxima  with  reference  to  the  electron 
interacting  with  molecular  surroundings.  We  have  not  yet  determined  that 
the  overall  system- -atom- atom  interactions  as  well  as  atom-electron 
interactions -- is  a  saddle  point.  It  is  entirely  possible  that  in  a 
reaction,  especially  in  a  reaction  in  solution,  the  electron  is  carried 
along  from  its  initial  to  final  state  in  a  local  energy  well  while  the  whole 
system  always  has  at  least  one  mechanically  unstable  degree  of  freedom.  The 
role  of  the  solvent  may  be  partially  to  solvate  the  migratory  electron.  The 
electron  in  the  transition  state  in  essence  would  be  trapped  between 
reactive  centers  and  molecules  of  solvent.  If  this  should  be  the  case  in 
some  instances,  some  of  what  has  been  learned  about  solvated  electrons  ought 


Co  apply  Co  Che  intramolecular ,  solvaced  eleccron  Cransfer  reaccion  [18]. 

On  Che  ocher  hand,  in  Chose  insCances  for  which  chere  is  a  clear  energy 
barrier  in  Che  CransiCion  sCaCe,  cunnelling  of  Che  eleccron  will  be 
necessary  Co  geC  ic  Co  Che  final  sCace.  Here  Chen,  Che  well -developed 
machinery  of  Che  diabacic  CransiCion  sCaCe  Cheory  should  apply  [1]. 

For  Che  pasc  30  years,  Che  eleccron  cransfer  Cheory  has  used  polar 
dieleccric  solvenC  model  as  an  inCegral  pare;  Che  oucer  sphere  solvenC  is 
usually  handled  in  varying  degrees  of  sophisCicaCion  and  complexiCy  as  a 
polar  conCinuum  [1].  Al Chough  local  molecular  vibraCions  are  frequently 
included  in  Che  analyses,  Chey  are  generally  examined  in  Che  harmonic  limit 
with  Che  machinery  of  Che  normal  mode  analysis.  The  overall  accuracy  of 
the  theoretical  accounC  of  observed  eleccron  transfer  is  a  matter  of  growing 
concern  [19] . 

It  is  possible  Chat  Che  dependence  of  Che  electron  Cransfer  reaccion  on 
local  vibrational  degrees  of  freedom  may  be  stronger  than  the  dependence  on 
Che  conCinuum  dielectric  response.  There  is  Che  well-known  example  of  the 
solvenC  independence  of  the  radiationless  transition  in  benzene  [20],  Energ 
loss  Co  vibrational  modes  is  known  to  involve  Che  C-H  stretch  motions;  there 
is  only  a  very  weak  dependence  on  solvenC.  Such  a  marked  insensitivity  Co 
Che  solvenC  may  not  apply  in  general  Co  electron  cransfer,  redox  reactions; 
ic  is  somewhat  early  Co  cell.  However,  Che  growing  difficulty  in  getting 
agreement  between  experiment  and  cheory  for  a  number  of  reactions  in 
several  solvents  poinCs  Co  a  more  important  involvement  of  local  vibrational 
modes.  The  role  of  Che  outer  modes  will  not  disappear,  of  course.  Their 
influence  can  be  incorporated  into  Che  local  vibrational  degrees  of  freedom 
as  mode  dressing  through  renormalization  [21] .  A  focus  of  attention  on  the 
local  modes  may  redress  Che  strong  emphasis  the  current  cheory  places  on  the 
outer  sphere  degrees  of  freedom.  The  strategy  for  determining  transition 


states  which  we  have  proposed  we  believe  is  well  suited  to  exploring  this 
issue . 

Our  purpose,  as  set  forth  in  this  note,  was  a  modest  one,  to  indicate 
that  it  is  possible  to  determine  transition  states  for  electron  transfer 
reactions  by  manipulating  the  center -of -gravity  of  the  mobile  electronic 
charge  density.  We  have  not  included  an  account  of  exchange  interactions  or 
configuration  interactions.  These  interactions  are  important,  and  will  be 
considered  in  due  course.  In  spite  of  the  implicit  computational  costliness 
of  the  approach  we  advocate,  we  believe  nevertheless  that  it  may  provide  a 
way  to  carry  out  useful  simulations  of  the  electron  transfer  reaction. 


This  work  was  supported  in  part  by  contracts  with  the  Office  of  Naval 
Research,  Arlington,  Virginia. 


Table  1 


Parameters  used  in  the  calculations: 
A.  Slater  basis  functions'*'10: 


-  1.0 


f.s  -  2.0 


Primitive  Gaussian  atomic  orbital'  : 

S  -  «  -  0.5  (Figs.  2-4):  S  -  0.4,  6  -  0.5  (Fig.  5) 

A  B  A  B 

Charges  and  interatomic  distances,  as  indicated  in  the  figures. 


Slater  functions  in  terms  of  a  scaled  Gauss  Ians  [16]  (t  -  term  symbol) 


Vr>  “  r  \ta.rr> 


Guassians  (cf.  ref.  16): 
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c.  Primitive  Gaussian  atomic  orbital: 
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1.  A  schematic  representation  of  the  charge  transfer  system.  The  atomic 
cores  are  labelled  A  and  B  and  are  located  on  the  z-axis.  The 

extent  of  the  underlying  valence  electronic  charge  is  indicated  by  the 
dotted  lines.  The  bold  dashed  line  indicates  the  extent  of  the 
transfer-electron  wavepacket. 

2.  The  energy  of  the  transfer  electron  in  the  field  of  two  neutral 
atoms  A:  the  reaction  is  A  +  A  •  A  +  A  ,  The  separation  between  the 
A- atoms  is  1.2  A. 

3.  The  same  system  as  in  Fig.  2.  The  interatomic  distance  is  2.5  A. 

4.  The  energy  of  the  electron  in  the  field  of  two  identical  cations: 

reaction;  A+  +  A  -  A  +  A+.  The  interatomic  separation  is  1.5  A.  The  height 
of  the  barrier  increases  with  increasing  A- A  separation. 

5.  The  energy  of  the  electron  in  the  field  of  two  dissimilar  cations  A  and 

B:  reaction;  A++  B  -  A  +  B+.  Here  6  -  0.4  A  and  S  -  0.5  A.  The 

•  b 


interatomic  separation  is  1.5  A. 
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